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Abstract 

We begin a systematic study of the enumerative combinatorics of 
mixed succession rules, which are succession rules such that, in the as- 
sociated generating tree, the nodes are allowed to produce their sons at 
several different levels according to different production rules. Here we 
deal with a specific case, namely that of two different production rules 
whose rule operators commute. In this situation, we are able to give 
a general formula expressing the sequence associated with the mixed 
succession rules in terms of the sequences associated with the compo- 
nent production rules. We end by providing some examples illustrating 
our approach. 

1 Introduction 

Among the many methods that have been developed to enumerate com- 
binatorial structures, the role of the ECO method has been growing in the 
last decade, thanks to its intrinsic simplicity and to the effectiveness of the 
combinatorial constructions it generates. The variety of problems in which 
the ECO method has shown its soundness ranges from enumerative and 
bijective combinatorics to random and exhaustive generation. 

The roots of the ECO method can be traced back to jCGHKj . where the 
authors study Baxter permutations and introduce for the first time the con- 
cept of a generating tree. Successively, West |WH IW2j introduced the notion 
of a succession rule to give a formal description of generating trees in the 
context of permutation enumeration and Barcucci, Del Lungo, Pergola and 
Pinzani [BDLPPJ extended the technique of generating trees, finding a gen- 
eral way of constructing combinatorial objects which can be often described 
using such formal tools. 

The classical ECO method (a detailed description of which can be found, 
for instance, in [BDLPP]) consists of a recursive construction for a class of 
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objects by means of an operator which performs a "local expansion" on the 
objects themselves. Typically, starting from an object of size n, an ECO 
construction allows to produce a set of new objects, of size n + 1, in such a 
way that, iterating the construction, all the objects of the class are obtained 
precisely once. If the construction is sufficiently regular, it can be often 
described by means of a succession rule, which is a system of the form 



The meaning of what is written above is the following. Each object 
of the class is given a label (k). When performing the ECO construc- 
tion, an object labelled (k) produces k new objects labelled, respectively, 
(ei (&:)), (e2(k)), . . . , (e^ik)). Moreover, the object of minimum size has label 
(a) (which is called the axiom of the succession rule). To have a graphi- 
cal description of a succession rule one usually draws its generating tree, 
that is the infinite, rooted, labelled tree whose root is labelled (a) (like 
the axiom) and such that each node labelled (fc) has k sons, labelled 
(ei(fe)), (e2(fc)), . . . , (ejfc(fe)), respectively. It is evident from this definition 
that we can introduce a notion of level on generating trees, by saying that 
the root lies at level 0, and a node lies at level n when its parent lies at level 



We remark that, from the above given definition, a node labelled (k) has 
precisely k sons. When a succession rule has this property it is often said to 
be consistent. However, one can also consider succession rules (and gener- 
ating trees) in which the value of a label does not necessarily represent the 
number of its sons, and this will be frequently done in the sequel. Moreover, 
we would like to warn the reader that, even if we will sometimes give defi- 
nitions using consistent succession rules (since this is the convention when 
working with the ECO method), we will constantly make use of succession 
rules which are not necessarily consistent. 

From the enumerative point of view, the main information encoded in 
a generating tree (and thus in its associated succession rule) is given by 
the level polynomial p n (x) = J2kPn,kX k , defined by setting p n ^ equal to the 
number of nodes labelled k at level n, and by the associated integer sequence 
(/ n )neNj which is defined, in terms of the level polynomials, as f n = p n (l), 
and represents the total number of nodes at level n. 

We point out that the infinite lower triangular array (p n ,fc)n,fceNi some- 
times called the AGT matrix [MV], or ECO matrix [FP2], often happens to 
be a Riordan array |Sp| . By definition, this means that every element p n & 
can be expressed by using a pair of formal power series (d(t), h(t)) in such a 
way that it is precisely the coefficient of t n of d(t)h(t) k . In this case, many 
counting properties of the generating tree can be found in an algebraic way, 
by using the related theory. 




(1) 



n — 1. 



2 



To give a quick example, consider the succession rule 



H -\ (k) -w (2) (3) (4) •■•(&)(& + 1) ' 

defining Catalan numbers 1, 2, 5, 14, 42, 132, . . . (see, for example, [BDLPPJ). 
The first levels of its generating tree can be depicted as follows: 




Here the level polynomials are p n (x) = Ylk^n,kX k , where b n ^ = 

2n-k+3 (^n+i" 3 ) are ^ e USU& 1 ballot numbers. Notice that, in this exam- 
ple, by shifting the column index k by two positions, it can be shown that 
B = On,fc)n,fceN = (C(t),tC(t)), where C{x) = is the generating 

function of Catalan numbers. Using the results of |Sp|, we can compute the 
row sums and the weighted row sums of B, thus obtaining the basic data for 
evaluating the probability distribution of the labels in the generating tree. 

The ECO method has been fruitfully applied to several problems, not 
only of an enumerative nature. For instance, using this technique it has been 
possible to develop efficient algorithms for the random [BDLP and exhaus- 
tive [BBGP] generation of combinatorial objects. Moreover, in |BGPP] the 
authors describe a general exhaustive generation algorithm (working for a 
wide family of structures) defining Gray codes not depending on the specific 
nature of the objects to be generated, but only on the properties of a succes- 
sion rule encoding an ECO construction for the objects under consideration. 
Due to its plentiful applications, it is then worth exploring in a deeper way 
the features and capabilities of such a method. 

Despite its wide range of applicability, there are many combinatorial 
constructions which cannot be naturally described by using the classical 
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ECO method (and classical succession rules) as exploited above. For in- 
stance, in |BM| a generalization of the method is considered, allowing suc- 
cession rules in which the labels are pairs of integers (rather than integers). 
Applications of this generalized method to the enumeration of pattern avoid- 
ing permutations are shown in the cited paper. 

A very strong limitation in the possibility of describing a combinatorial 
construction by means of a succession rule lies in the definition we have 
given of this fundamental tool. The generating tree of a succession rule has 
the property that, if the level of a node is n, then the level of all its sons 
is n + 1. From a combinatorial point of view, this means that a (classical) 
ECO construction performed on an object of a given size produces objects 
of the successive size. However, it may well happen that a combinatorial 
construction, having all the reasonable features to be called ECO, does not 
behave in the standard way with respect to the notion of size. More precisely, 
starting from an object of size n, we can construct new objects whose sizes 
are greater than n (but not necessarily equal to n + 1). The formalization of 
these concepts leads to the notion of what can be called a mixed succession 
rule. Roughly speaking, the idea is to consider a set of (possibly different) 
succession rules acting on the objects of a class and producing sons at differ- 
ent levels. To be more formal, we introduce here the simplest instance of this 
general situation, by considering two succession rules producing their sons 
at the two successive levels. These will be called doubled mixed succession 
rules. Given two succession rules £1 as in (pQ) and 



we define the doubled mixed succession rule associated with the pair (f2, S) 
with axiom (c) to be the succession rule (c)f2 +1 E +2 , defined by 



The generating tree associated with (c)f2 +1 £ +2 has the property that 
each node labelled (k) lying at level n produces two sets of sons, the first 
set being (ei(fc)), (e2(fe)), . . . , (ek(k)) at level n + 1 and the second one being 
(di(k)), (d2(k)), . . . , (d k (k)) at level n + 2 (so that it produces a total of 2k 
sons). 

To justify our interest in this kind of notion, we remark that instances of 
(general) mixed succession rules have occasionally surfaced in some previous 
works; to cite only one example, in |GPPj vexillary involutions are enumer- 
ated by making use of a specific mixed succession rule. The first systematic 
treatment of mixed succession rules has been undertaken in |FPPR2j , where 





(c)ft +1 £+ 2 : { (k) 



+1 



( ei (k))(e 2 (k))-..(e k (k)) 
(d l (k)){d 2 {k))---{d k {k)) 



+2 
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the special case £ = $7 has been examined in great detail. The present paper 
represents the first attempt to tackle the general case, aiming at developing 
a general theory of mixed succession rules. More precisely, the main goal 
would be to succeed in expressing the sequence associated with a mixed suc- 
cession rule (c)0 +1 S +2 in terms of the sequences associated with Q and S, 
possibly changing the axioms. For this reason, in section [3] we study some 
enumerative properties of what we have called production rules, which are, 
by definition, succession rules without the axiom. The problem of studying 
mixed succession rules in its full generality seems quite difficult; it is some- 
how related to the theory of power series in several noncommuting variables. 
Here we deal only with a special case, namely when the two rule operators 
(see section [2]) of and S commute. In this situation, we are able to find a 
general formula for the sequence associated with (c)f2 +1 £ +2 ; moreover, we 
also describe some examples our theory can be applied to. 

We would like to remark that this problem has also been considered from 
the point of view of Riordan arrays [BMSJ . Each Riordan array determines 
a specific sequence (a^gN, called the A-sequence of the array, such that, 
for every n, k G N: 

oo 

b n +i,k+i = ^] a-kb n ,n+k- 

k=0 

When the ^-sequence contains integer numbers only, it is related to the 
succession rule of an associated generating tree (if any) as shown in [MVj. 
However, it can happen that the A-sequence has a complicated expression, 
whereas the ^-matrix (as defined in |BMSj ) is simple. This corresponds to 
an ECO construction in which the elements of size n produce objects of 
different sizes (greater than n + 1). 

In closing this introduction, we recall some notations we will frequently 
use in the next pages. 

The sets of natural and real numbers will be denoted N and R, respec- 
tively. 

The following linear operators on the vector space of one-variable poly- 
nomials will be often considered: x (respectively, t) is the operator of multi- 
plication by x (respectively, t), D is the usual derivative operator, and T is 
the factorial derivative operator, which is, by definition, the linear operator 
mapping x n into 1 + x + ■ ■ ■ + x n_1 = ^7=0 x% (f° r n — 1) anc ^ 1 0- 

2 Preliminaries on rule operators 

Given a succession rule £1 as in (pQ), we can associate with it a linear 
operator on the vector space of one- variable polynomials to be denoted 
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L = Lq (the subscript will be omitted when it is clear from the context). 
To define such an operator, we use the canonical basis (x n ) ne -N'- 

L : K[x] — ► K[x] 
1 i — ► x a 

x k i — ► x ei{k) H h x 6fe(fc) , if k appears in Q,, 

x h i — ► hx h , otherwise. 

As it is easy to understand, each enumerative property of $7 can be 
suitably translated into some property of L. For instance, if (/ n )neN is the 
numerical sequence associated with fi, for any n £ N the following equality 
holds: 

f n =[L n +\l)] x=1 , 

where we have used square brackets to denote the operator of evaluation at 
a specific value. The linear operator L is called the rule operator associated 
with tt. We refer the reader to IFPPRIIIFTT^IFTTIIFPI] for the definition, 
the main properties and some applications of this notion. 

Remarks. 

1. In the above definition we have denoted with K a generic fields of 
coefficients. For what concerns us here, it is largely enough to take K = 
R. Actually, all the theory of rule operators could be equally developed 
on the semiring module N[x] of polynomials with nonnegative integer 
coefficients. 

2. If p n (x) = ^2kPn,k xk is the n-th level polynomial of it is 

Ln(p n (x)) = p n+ i(x). 



Examples. Here are some examples of rule operators associated with 
more or less well-known succession rules. All these examples can be found, 
for example, in |FPlj . 

1. A succession rule for factorial numbers: 

\ (k) * (k + i) k ' 

associated rule operator: 

L(l) = x 
L(x k ) = (x. 2 D){x k ) = kx k+ \ k>l. 
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2. A succession rule for arrangements: 

~ (2) 



(fc) (k)(k + l) k - 1 ' 
associated rule operator: 
L(l) = x 2 

L(x fe ) = (x 2 D-x+l)(a; fc ) = x fc + (fc - l)x fc+1 , fc > 1. 
3. A succession rule for involutions: 



(fc) -w (fc — l) fc_1 (fc + 1) ' 
associated rule operator: 
L(l) = x 

L(x fc ) = (D + x-x- 1 )(x fc ) = (ife-l)x fc - 1 +x fc+1 , fc>l. 

4. A succession rule for Bell numbers: 

(1) 



(fc) ~» (fc) fc ~ 1 (fc + l) ' 
associated rule operator: 
L(l) = x 

L(x fc ) = (x£> + x- l)(x fc ) = (fc- l)x fe + x fe+1 , fc>l. (2) 

Further examples involving the factorial derivative operator T can be 
considered. 

5. A succession rule for Catalan numbers: 

(1) 

(fc)-(2)(3)---(fc)(fc + l) ' 
associated rule operator: 

L(l) = x 

k+i 

L(x k ) = (x 2 T)(x fc ) = ^2x\ fc>l. (3) 

i=2 
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6. A succession rule for Motzkin numbers: 
(1) 

(1) - (2) ; 
(k) -w (1)(2) •••(& — 1)(A; + 1) jfe>l 

associated rule operator: 

L(l) = x 

k-i 

L{x k ) = (xT + x- l)(x fc ) = ^V + x fc+1 , fc>2. (4) 

i=l 

3 Production rules 

With the expression production rule we will mean here a succession rule 
without its axiom. Hence the generic form of a production rule is 

(k)^( ei (k))(e 2 (k))...(e k (k)). (5) 

Clearly, in the same way a succession rule determines a unique numer- 
ical sequence, a production rule defines a family of sequences (/n )neN> 
depending on the axiom (a) which we choose for the rule ((5]). 

From now on, given a succession rule O as in ([T]), we will denote with L a 
the associated rule operator, (a) being the axiom of f2. Using this terminol- 
ogy, given a production rule as in ([5]), the family of operators (L a ) ag N will 
be called the family of rule operators associated with the production rule. 
In this section we will be interested in finding formulas to relate the various 
sequences associated with the same production rule. To this aim, we start 
by observing the following (very easy but quite important) facts. 

1. For any a, b 6 N, we have L a (x k ) = Lf,(x k ), when k ^ 0, and L a (l) = 
x a , Lfc(l) = x b . For this reason, in what follows, for k ^ 0, we will 
simply write L(x k ), without specifying the axiom, and we will speak 
of the rule operator associated with the production rule whenever we 
restrict our attention to the subspace X-ftTfx] (i.e. the subspace spanned 
by the positive powers of x). 

2. The n-th term of the numerical sequence of the family with axiom (6), 
that is fn \ can be computed using the following formula: 

/« = [L n (x b )] x=1 . 

In the same way, to compute fn +l ^ we get: 

= [L n (x b+1 )] x= i = [L n x(x b )] x=1 . 
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Then, if one knows the Pincherle derivative of L, which is, by defini- 
tion, the operator L' = Lx — xL (see |RKOj ). it should be possible 
to express the operator L n x as a linear combination of monomials of 
the kind x a L@. This should allow, at least in principle, to obtain an 
expression for/i, 6+1 ^ in terms of known quantities (namely fm +1 \ with 
m < n, and , with a < b). To better understand how to proceed in 
the concrete cases, it is convenient to have a look at a specific example. 

3.1 A Bell-like production rule 

Let us consider the following production rule: 



On the nonzero powers of x, the rule operator associated with u acts as 
follows: 

L(x k ) = (xD + x- l)(x k ), Vfc^O. 

Thus, on we have that L = xD + x— 1. Concerning the Pincherle 

derivative L' , on xET[x] we have: 

Lx = (xD + x - l)x = x(xD + 1) + x 2 - x = x 2 D + x 2 , 
xL = x(xD + x - 1) = x 2 D + x 2 - x, 

(here we have used the well-known identity D' = Dx — xD = 1). Therefore, 
we get: 



Our next goal is to express the operator L n x in such a way that some 
identity between the terms of the sequences associated with u can be deter- 
mined. 

Proposition 3.1 For any n > 1, we have: 



Proof. First observe that, from ([6]), we get Lx = x(l + L), whence: 



(k) ~» (jfe^-^ib + l). 



(6) 




k=0 




x(l + L) 



n 



x(l + L)(l + L) 



n-1 



Lx(l + L) 



n-1 



and then, by iterating: 



x(l + L) n = L n x 
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The above formula is easily seen to hold for n > 0. By iterating, we then 

get: 

L n x = x(l + L) n = x(l + L) n - 1 +xL(l + L) n - 1 
n-l 

-k 



fc=0 v 7 

u—n \ / 



Formula ([7]) yields an interesting recursive expression for the family of 
sequences associated with L. 

Corollary 3.1 Denoting by (fn^)neN the sequence associated with L with 
axiom (b), it is 



f (6+1) _ ST f n ~ A f W , f 

Jn - 2^1 \ U Jk+l^Jn-1 ■ 
1,—n \ / 



k=0 

Proof. Indeed, using the above proposition, we get: 

n-l 



/ — 1 \ 

= [L n x(x b )} x= i = J2[ k )[^L k+1 (x b )] x=1 + [L^(x b )} x=1 



k=0 

as desired. ■ 

As a further consequence, one can, for example, determine the generat- 
ing function f^ b \x) of (/ji )neN- Standard generating function arguments 
provide the following result. 

Corollary 3.2 Denoting by B the binomial transform operator on exponen- 
tial generating functions (such as, given an exponential generating function 
f(x), it is B{f{x)) = e x f{x)), we have: 

f^(x) = B b - 1 (fW(x)). 

Since f^'(x) = e e *~ l is the exponential generating function of Bell numbers 
1,1,2,5,15,52,203,..., we get 

/W( x ) =e «-+(6-D*-l 
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3.2 Further examples 

It is obvious that what we have done for the rule operator xD + x — 1 
can be done (at least in principle) for any other rule operator. In particular, 
we mention here two further cases, giving only the final results, and leaving 
to the reader the details of the proofs. 

1) A Catalan-like production rule. Consider the production rule 

u: (fc)~»(2)(3)---(fc)(fc + l), 

which is related to Catalan numbers. Indeed, it is known that the 
generating function of the sequence with axiom (b) is C b (x) (recall 
that C(x) is the generating function of Catalan numbers). We can use 
our approach to rediscover this result. 

Indeed, we first observe that L = x 2 T is the rule operator associated 
with to (as usual, we restrict our attention to the subspace xK[x]). We 
start by computing the Pincherle derivative L'. 

Lemma 3.1 Let 02 : xK[x] — ► x_RT[a;] be the linear operator defined 
on the canonical basis by setting C2(x n ) = x 2 , for every n > 0. Then 
we have that 

L' = c 2 . 

From this lemma we obtain the following result, which is essential in 
deriving the final recursion. 

Proposition 3.2 For every n G N, it holds: 

n-1 

L n x = *L n + Y J L i c 2 L n - 1 - i . 

i=0 

Finally, after having observed that, for all i G N, it is clearly 

(2) 

[L l C2(p(x))] x= i = fl ■ p(l), we get the following recursion for the 
family of sequences associated with the starting production rule. 

Proposition 3.3 For every b, n G N, we have: 

n-1 

f (6+1) = f(b) , V f (2) Ab) 

Jn Jn ' / j Jj Jn—l—i' 

As a consequence, the generating function of (fn^)ne'N is given by 
the convolution C b (x) of the generating function C(x) = 1 ~^~ 4x of 
Catalan numbers. 
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2) A Motzkin-like production rule. Consider the production rule 

u: (A0~»(l)(2)---(fc-l)(fc + l), 

which is related to Motzkin numbers. Namely, it is known that u, with 
axiom (1), yields the sequence of Motzkin numbers 1, 1, 2, 4, 9, 21, 
The rule operator associated with oj is L = xT + x — 1 (also in this 
case we refer to xK[x]). We can proceed analogously as we have done 
in the previous example. 



Lemma 3.2 Let c\ : xK[x\ — ► iJfJi] be the linear operator defined 
on the canonical basis by setting ci(x n ) = x, for every n > 0. Then we 
have that 

U = ci. 



Proposition 3.4 For every n GN, it holds: 

n-1 

L n x = xL n + ^L i ciL"- 

i=0 



Proposition 3.5 For every b, n £ N, we have: 

n-1 



f (6+1) = f(b) , V f (1) f {b > 
Jn Jn > / j J i J n—1- 



y(i) m 

i=0 

As a consequence, the generating function of (/n^) n eN is given by 
the convolution / (b) (x) = f {h ^ 1 \x)(xM{x) + 1), where M{x) = 
1 ~ x ~^2x 22X ~ 3x2 ^ s the generating function of Motzkin numbers, and 



so f( b )( x ) = M(x) ■ (xM(x) + 1) 



6-1 



4 Commuting rule operators 

In this section we finally tackle the problem of finding an explicit expres- 
sion for the numerical sequence associated with a mixed succession rule. As 
we have said in the introduction, in what follows we will deal with the sim- 
plest case of a mixed succession rule, namely the case of a doubled rule: this 
means that, in the associated generating tree, each node at level n produces 
a set of sons at level n + 1, according to a production rule u>, and another 
set of sons at level n + 2, according to another production rule a. If (6) is 
the axiom of the doubled rule, such a generating tree can be synthetically 
represented as follows: 

As we have declared in the introduction, our main aim is to determine 
an expression for the sequence associated with the doubled rule (with axiom 
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(b) 



level n 



level n+1 
level n+2 



(b)) in terms of the sequences associated with the production rules u and 
a. Let's start by fixing some notations. First of all, L and M will be the 
rule operators associated with u and a, respectively. Using production rules, 
our doubled mixed succession rule will be denoted (b)co +1 a +2 , whereas in 
terms of rule operators it will be {b)L +l M +2 . Moreover, we will denote p n {x) 
the level polynomials of {b)uj +l a +2 . Finally, f^ b \x,t) will be the bivariate 
generating function of the generating tree, where t keeps track of the level 
and x keeps track of the label. Our first result is an expression for f®(x,t) 
in terms of the rule operators L and M. 

Proposition 4.1 Denoting by ~ 1 the compositional inverse of an operator, 
we have: 

f {b) (x,t) = (l-tL-t 2 M)- 1 (x b ). 

Proof. The argument to be used here is analogous to the one used in 
[FPPR2| for jumping succession rules. Since each node at level n can be 
generated either by a node at level n — 1 (according to to) or by a node at 
level n — 2 (according to a), we have the following expression for p n (x): 

p n {x) = L(p n _i(x)) + M(p n _ 2 (x)). 

If we impose, by convention, that Pi(x) = 0, for i < 0, then the above 
expression is meaningful when n > 1 (recall that, under our assumptions, 
Po(x) = x b ). In order to translate the above recursion into generating func- 
tions, we multiply by t n both sides of the above equality and sum up for 
n > 1, thus obtaining: 

Y^Pn{x)t n = L(p n ^(x))t n + M(p n _ 2 (x))f\ 

n>l n>l n>l 
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whence, using linearity: 

n>l \n>l ) \n>l ) 

Since f( b \x,t) = ^2 n >o Pn(x)t n , we will then get: 

/(*>(*,*) - x fc = tL(/( b )(x,t)) +t 2 M(f<- b \x,t)), 

whence 

/^(x,*) = (l-tL-t 2 M)- 1 (s b ). ■ 
Expressing the operator (1 — tL — t 2 M) _1 using power series, we have 

that 

fW(x,t) = ^2t n {L + tM) n (x b ). (8) 

n>0 

Therefore, it is now clear that, if we want to know the sequence associ- 
ated with the doubled rule, we need to find an expression for the binomial 
(L + tM) n . In general, this is a nontrivial problem, since the linear opera- 
tors L and M usually do not commute. We are then led to first take into 
consideration just a special class of pairs of rule operators. More precisely, 
in the rest of the paper, we will assume the following hypothesis: 

L ed M commutes, i.e. LM = ML. 

Using an algebraic terminology, it is said that the commutator [L, M] = 
LM — ML is equal to zero. 

Now let's come back to our problem, that is the determination of an 
expression for the binomial (L + tM) n . We have the following, crucial result. 

Theorem 4.1 Denoting by (Jf\x) = ^iti'fiX' 1 the r-th level polynomial 
of the generating tree of a with axiom (s) and by (li^)ngN the numerical 
sequence associated with to with axiom (a), if (/n^) n eN is the sequence de- 
termined by (b)uo +1 a +2 , we have: 



k>0 V 7 i 



(9) 
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Proof. From the fact that [L, M] = we get immediately: 

(L + tM) n = (?) t k L n ~ k M k . 
k=o ^ ' 

As a consequence, equality © can be rewritten as: 
f®(x,t) = ^t n ^t fc ^L n - fc M fc (x 6 ) 

n>0 k=0 ^ ' 



J2\J2i n I k )L n - 2k M k (x b ; J t". 



n>0 \fc>0 



From the above expression we immediately deduce that p n (x) = 
Efc>o { n ~ k k )L n - 2k M k {x b ), and so the n-th term of the sequence associated 
with the doubled rule, which is fn = Pn(l)j can be computed as follows: 



p„(i) 



E 

fc>0 



n — k 
k 



L n ~ 2k M k (x b ) 



k>0 



x=l 



x=l 



Ef": fc )E^' 



fc>0 



(&),« 
n-2fc' 



and this is precisely our thesis. ■ 

Therefore we have succeeded in finding a formula expressing the nu- 
merical sequence associated with a doubled mixed succession rule when the 
related rule operators commute. Specifically, our formula involves: 

• the distribution of the labels of the production rule a with axiom (b) 
inside its generating tree (i.e. the coefficients ^\)\ 

• the sequences associated with the production rule to (i.e. the coeffi- 



cients I 



n-1k>- 



5 Examples 

We close by giving two applications of formula Q. The first case is 
somehow trivial (but leads to interesting enumerative results) , since we deal 
with the identity operator, which does not raise any problem concerning 
commutativity. However, in any other case, we need to determine some pairs 
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of commuting rule operators. To this aim, the easiest way is perhaps to 
fix a rule operator L and then find the general form of the rule operators 
commuting with L, which is precisely what we have done in our second 
example. 



5.1 The identity operator 

As it is obvious, the identity operator 1 commutes with any linear op- 
erator. Therefore, if L, M are any rule operators, we can consider the two 
doubled mixed succession rules (6)L +1 1 +2 and (6)1 +1 M +2 . Let us analyze 
the two cases separately. 

Consider first (6)L +1 1 +2 . To apply theorem 14. 1\ we observe that 

(^)neN is the sequence determined by L with axiom (s), whereas {%) is 
the r-th level polynomial of the succession rule determined by the identity 
operator 1 with axiom (s), and so it is trivially iif\x) = x s , whence 

(a) f 1 , i = S 

*V 10, i ^ s ' 



Thus, denoting by (/f ) n eN the sequence determined by the doubled 
mixed rule (6)L +1 1 +2 , we get: 



Jn — I ^ ) L n-2k- 

k>0 



Moreover, if f( b \x) and l^ b \x) are the two generating functions of the 
above sequences, standard arguments leads to the following result. 

Corollary 5.1 For the generating function f^ b \x) we have: 
Examples. 

1. If L is the rule operator of Catalan numbers described in ([3]), the 
sequence determined by the mixed rule (1)L +1 1 +2 is sequence A105864 
in [Si], which has no significant combinatorial interpretation. In order 
to find one, we consider a special class of parallelogram polyominoes. 
A 1-2 column parallelogram polyomino is a parallelogram polyomino 
whose cells can be either monominoes or dominoes, such that 

i) each column is entirely made either of monominoes or of domi- 
noes; 
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ii) given a set of consecutive columns starting at the same height, 
the leftmost one must be made of monominoes. 

Such a class of polyominoes can be constructed as follows, according 
to the semilength. Given a polyomino V of semilength n and such that 
its rightmost column has k — 1 cells, we construct the following set of 
new polyominoes: 

• either add a new rightmost column made of monominoes ending 
at the same height of the rightmost column of V, or 

• add a new cell on the top of the rightmost column of V (such a 
new cell will be a monomino or a domino according to the type 
of the column), or 

• add a new rightmost column made of dominoes, starting and 
ending at the same heights as the rightmost column of V . 

An instance of this ECO construction is shown in figure [TJ 



+ 1 

3» 



+ 1 

5* 



+ 1 

- 



+ 1 

5* 



Figure 1: Our ECO construction performed on a polyomino of semilength 
14 and such that the rightmost column has 3 cells. 

As it is clear, in the first two cases polyominoes of semilength n + 1 are 
produced, whereas in the third polyomino of semilength n + 2 

comes out. Now the reader can check that the above construction can 
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be encoded by the following succession rule: 



(1) 

(k) 



3(2)(3)...(fc)(fc + l) , 



if (k) 



that is precisely the mixed succession rule (1)L +1 1 +2 . 

2. Taking for L the rule operator of Motzkin numbers recalled in (J3J) 
(and choosing again (1) as axiom), we get sequence A128720 of [Si] . 
One of the given combinatorial interpretations for such a sequence is 
the following: it counts the number of 2- generalized Motzkin paths, i.e. 
paths in the first quadrant from (0,0) to (n, 0) using steps U = (1, 1), 
D = (1, —1), h = (1,0), and H = (2,0). Various kinds of generalized 
Motzkin paths have been extensively studied in the literature, see for 
example |dMSl ISu| . The mixed succession rule arising in this case is 
the following: 



(1) 

(*) 



(l)(2)...(fc-l)(fe + l) 



if (k) 




Figure 2: The generating tree of Q 



It is interesting to notice that f2 (whose generating tree is depicted in 
figure Ej) indeed describes an ECO construction for the above class of 
paths. Leaving the details to the interested reader, we quickly justify 
this claim: given a 2-generalized Motzkin path, consider its last de- 
scent, i.e. the final sequence of the path free of U steps. Construct a 
set of new paths as follows: either replace each h step with a U step 
and add a D step at the end, or just add a h step at the end, or sim- 
ply add a H step at the end. In the first two cases the length of the 
path is increased by 1, whereas in the last case it is increased by 2. 
It is now easy to show that this construction is encoded precisely by 
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the mixed succession rule £1. We also remark that another interpreta- 
tion of sequence A128720 is provided in [SI], namely using skew Dyck 
path. It would be interesting to use the above mixed rule to describe 
a construction for this latter combinatorial structure as well. 

To conclude this example, we also notice that, according to [BMSj, the 
rule f2 is alternatively described by the ^-matrix 

1 •• 

1 1 1 1 •• 

which implies that the associated ECO matrix is actually a Riordan 
array. Its ^4-sequence has generating function: 



1 - t + t 2 + Vl - 2t + 7t 2 - 10t 3 + 5t ? 
A[t) ~ 2(1— t) 

= 1 + 2t 2 + t 3 - t A + 6i 6 + 5t 7 - 16t 8 + ■■■ , 

and this shows that a direct dependence of row n + 1 from row n is 
very unlikely. Using the theory of Riordan arrays, we can determine 
the formal power series d(t) and h(t) defining our ECO matrix. More 
precisely, since h(t) = tA(h(t)), we find: 



h(t) 




l-3t-t 2 \ 
1 + t-t 2 J 

= t + 2t 3 + t 4 + 7t 5 + 10t 6 + 37t 7 + 82t 8 + ■■■ . 

Since column is not privileged, we have d(t) = so our Riordan 
array is completely determined. Denoting by S(t) = ^2 n >oS n t n the 

generating function of the row sums of the array, since S(t) = j^rfhi , 
we get: 



_ i-t-t 2 -^{i-t-t 2 )(i-?,t-W) 

b[t} - 2T 2 

= 1 + t + 3t 2 + 6t 3 + I6t 4 + 40t 5 + 109t 6 + 297t 8 + 



Using standard methods of asymptotic analysis, we obtain the asymp- 
totic value: 

K ip- n ~ 2 /2n + 4 N 



iP~ n - 2 (2n + 4\ 
4"+ 2 (2n + 3) \ n + 2 )' 



where K = ^22^13- 78 « 1.149838276... and ip = ^|=^. For n = 60 
it is Set) ~ 4.960467337... x 10 28 , whereas the above formula gives 
4.949459297... x 10 28 , with a relative error of 0.22%. 
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3. 



If L is as in ([2]), defining Bell numbers, the resulting sequence (/n ) n eN 
starts 1,1,3,7,22,75,... and is not recorded in [Si] . Thanks to our 
theory, it is possible to give a combinatorial interpretation to such a 
sequence, by performing an ECO construction described by the mixed 
succession rule: 



n : < 



(1) 

(*) 



+1 

+2 



(*)*" 
(*) 



We call a lacunary partition of [n] = {1,2, ...,n} any partition of 
a subset S of [n] such that [n] \ S is disjoint union of intervals of 
even cardinality. For instance, the partition {{1, 8, 12}, {2}, {3, 9}} is 
a lacunary partition of [14]. An ECO construction for the class of 
lacunary partitions works as follows: given a lacunary partition ir of 
[n], construct a set of new lacunary partitions by either adding the 
block {n + 1}, or adding n + 1 to each of the block of it, or else leaving 
7r unchanged, but thinking of it as a lacunary partition of [n + 2]. Of 
course, performing one of the first two operations leads to a lacunary 
partition of [n + 1], whereas the last one produces a lacunary partition 
of [n+2] . The reader can now check that such a construction is encoded 
by Q. 

Now consider (6)1 +1 M +2 . In this case, the sequence (Zn )nsN is the one 

(s) 

determined by the identity operator, and so l n = 1, for all n S N. On 
the other hand, the polynomial fii (x) is the level polynomial of the rule 
associated with M with axiom (s). Applying theorem I4.1|. for the sequence 
) nS N determined by (6)l +1 Af +2 we then get: 

fc>0 ^ ' i k>0 V ' 

where, of course, [m n ) n &i is the sequence associated with the rule operator 
M with axiom (s). 

As in the preceding case, we can easily recover the generating function 
of the sequence f n b ^ starting from that of the sequence rrin ■ 

Corollary 5.2 For the generating function f^ b \x) we have: 

Examples. We leave as an open problem that of finding an ECO con- 
struction described by (1)1 +1 M +2 for each of the structures mentioned in 
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the next examples. In the first of them, we also provide a rigorous proof 
(via Riordan arrays) of the fact that a certain sequence on [SI] comes out, 
whereas the details of the second example are left to the reader. 

1. If M is the rule operator of Catalan numbers described in (|3|), the se- 
quence determined by the mixed rule (1)1 +1 M +2 is sequence A090344 
in [Si] . Such a sequence counts Motzkin paths without horizontal steps 
at odd height. 

From an enumerative point of view, observe that (1)1 +1 M +2 can be 
explicitly written as: 



(1) 

(k) 



+1 

+2 



(k) 



(2)(3)--- (*)(* + 1) 
According to [BMSj . this corresponds to the A-matrix 



11111 
10 



which gives rise to the following vertically stretched Riordan array (see 
[CMS] ): 



n\k 


1 


2 


3 


4 5 





1 








1 


1 








2 


1 


1 






3 


1 


2 






4 


1 


4 


1 




5 


1 


7 


3 




6 


1 


13 


8 


1 


7 


1 


24 


18 


4 


8 


1 


47 


40 


13 1 



The entries m n ^ of the above array obey the recursion m n+ \^+i = 
m n,k+i + Y^jLo m n-l,k+j- To obtain a proper Riordan array (to be 
denoted (d(t), h(t))), we can simply shift column fe+1 up by A; position 
(for all k), so that the entries of the new (lower triangular) array 
will now satisfy the equalities p n +i,k+i = Pn,k+i+Y,'jLoPn-j,k+j- These 
recurrence relations are described by the (infinite) ^4-matrix: 



/ 



\ 



1 

1 

10 

10 

\ i i o o o 
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To get the generating function of the ^4-sequence of our Riordan array 
we can use the equality A(t) = t j A(t)-ipW(t) (shown in |BMSj ). 

where P^(t) is the generating function of row j in the ^4-matrix. In 
our case, we have Pt°](t) = 1 + t and pW(t) = V , for j > 0, which 
yields: 



A(t) 



1 + 1 + tA(t)~H + t 2 A(t)~ 2 t 2 + t 3 A(t) _3 t 3 + 



t + l + 



+ 



+ 



+ 



A{t) A(t) 2 A(t) 3 
The solution of the above equation is: 

1 + t + t 2 + VI + 2t + 3t 2 - 2t 3 + t 4 



i + 



A-t 2 ' 



A(t) 



1 + i + t 2 - t 3 + t 4 - 2t 6 + 4i 7 - 3t 8 - 5t 9 + • • • . 



Now we are ready to compute d(t) and h(t). Using the well known 
formula h(t) = tA(h(t)), we get an equation of degree two, whose 
solution is: 



h(t) 




l-t- At 2 



l-t I 

t + t A + 2t° + 3t 4 + 6i 5 + lit 6 + 23t 7 + A7t 8 + 102t 9 + 



and obviously d(t) = j—r. Our original (stretched) Riordan array is 
therefore: 



1 



1 



l-t' 2 



1 



1 - t - At 2 



l-t II 

and the generating function of the sequence of its row sums is: 

oo 



S(t) = ^2 S n t n = - 



d(t) 



1 



n=0 



h(t) 2t 2 



1 



l-t- At 2 

T^t 



It is possible to find an asymptotic value for S n . The equation 1 — 
t - At 2 = has the two solutions r% = ~ v/ ^ +1 -0.6403882032... 

and r 2 = ^±J^I » 0.3903882032..., and therefore this last value is the 
dominating singularity. Using Bender's theorem |BJ, we obtain: 



1 



l-t- At 2 

T^t 



l 

~2 
K 



i - tjn 
l-t 

l 



[t n+2 ]Ji 



- t=T2- 

2n + A 



t 

r-2 



2 (2n + 3)r% +2 A n + 2 \n + 2 
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In the above formula, we have used the notation [t n ]f(t) to denote 
the coefficient of t n in the formal power series fit). For instance, for 
n = 50 we have S50 ~ 1.091877333... x 10 18 , while the approximate 
value is 1.075272279... x 10 18 , with a relative error of 1.54%. Since 
r\ < 1, its contribution decreases as n increases, but for small values 
of n it cannot be ignored. 

2. Taking for M the rule operator of Motzkin numbers recalled in (j4j) 
(and choosing again (1) as axiom), we get sequence A026418 of [SI]. It 
counts ordered trees having no branches of length 1, according to the 
number of edges. 

3. If M is as in ([2]), defining Bell numbers, the resulting sequence 
(/PjneN starts 1,1,2,3,6,11,23,47,103,... and is not recorded in 



5.2 A factorial-like rule operator 

Consider the rule operator L = x 2 L> associated with the production rule 



We start by determining the family of sequences related to L. 

Lemma 5.1 If (4?) n gN is the sequence determined by L with axiom (b), 
then we have, for all n 6 N; 



where (x) y = rr(x — 1) • . . . • (x — y + 1) denotes the usual falling factorial. 

Proof (sketch). Use a simple induction argument. For b = 1 it is well- 
known (see, for instance, |FP1] ) that In = n - Now observe that the recur- 
sion defined by the associated production rule implies that = bln +1 \ 
whence it is easy to derive the thesis. ■ 

According to our program, we start by computing the general form of a 
rule operator commuting with L. 

Theorem 5.1 Let M be a rule operator such that M(l) = x a , for some 
a G N. Then, M commutes with L if and only if 



[SI]. 



(k) ~* (k) k+1 . 




M = L [a] 



(a -1)! 



£> a x' 
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Proof. Suppose that M commutes with L. On the polynomial 1 it is 

M{x) = M(x 2 L>(l)) = x 2 D(M(l)) = x 2 D(x a ) = ax a+1 . 

Now suppose by induction that M(x n ) = n( a '^7 )x a+n . Then, on x n+l 
we have: 

M(x n+1 ) = -M(nx n+1 ) = -M(x 2 D(/)) = -x. 2 D(M(x n )) 
n n n 

We have thus showed that M(x n ) = n{ a+ a n _~ 1 )x a+n , that is M = L [a] = 
v +1 r D a x a ~ 1 , as desired. 



(^T)T 

As far as the converse is concerned, we leave to the reader the proof of 
the fact that the operator L[ tt j = (a-i)! ^° x " 1 commutes with L. ■ 

Consider now the case a = 2, so to obtain the operator M = Lr 2 i = 
x 3 Dx 2 . The mixed succession rule Qb = {b)L +1 M +2 is 

' (6) 

Q b : < (k) il (k + l) k . (10) 

+| (jfe + 2)*(*+i) 

The first levels of the corresponding generating tree, when 6 = 1, appear 
as in figure 3. 




Figure 3: The first levels of the generating tree of (1)L +1 M +2 . 

In order to apply theorem 14.11 we need to know the sequence iff deter- 
mined by the production rule associated with L with axiom (s) and the level 
polynomials /J*i (x) of the generating tree related to M with axiom (s). 

The first information is provided by lemma [5TTI that is iff = n '( n ^^ 1 )- 
As far as the polynomials (x) are concerned, we observe that, in the 
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generating tree associated with M, only one label appears at any given level; 
more precisely, the only label at level r is (s + 2r). Therefore iif\x) consists 
of only one monomial, and we simply have to determine its coefficient. The 
following, simple lemma finds this coefficient. 

Lemma 5.2 The generating tree associated with M having axiom (s) has 
(s) 2r nodes at level r (each of which is labelled (s + 2r)), where (x) y = 
x(x + f) • ... • (x + y — f) denotes the raising factorial. 

Proof. At level and 1 there are, respectively, 1 and s(s + l) nodes. By 
induction, suppose that at level r we have (s) 2r nodes; since each of them 
is labelled (s + 2r), it produces (s + 2r)(s + 2r + 1) sons, whence the thesis 
immediately follows. ■ 

As a consequence, we have that fir (x) = (s) 2r x s+2r , which means that 

/4*i+2r = ( s ) 2r > whereas = 0, for j / s + 2i. 

We are now ready to apply theorem 14.11 thus getting for the sequence 
(/ft neN determined by (b) L +1 M +2 the following formula: 



where (-F n )neN is the sequence of Fibonacci numbers. 

To give a combinatorial interpretation for the sequence (/n ) n eN we 
refer to the generating tree of the mixed rule fife. In what follows, S n will 
denote the set of coloured permutations of [n], i.e. permutations whose ele- 
ments can possibly be coloured (a coloured element will simply be overlined). 
Moreover, we introduce the notion of paired coloured permutation, to mean 
a coloured permutation such that, if / denotes the set of coloured elements, 
then / is a disjoint union of intervals of even cardinality. So, for instance, 
the permutation 931465728 is a paired coloured permutation belonging to 



Proposition 5.1 Given b G N, fix a permutation r G S&_i. Then /„ is the 
number of paired coloured permutations it G Sb +n _\ in which the elements 
1, 2, . . . , b — 1 are not coloured and appear in tt as a pattern isomorphic to 



Proof. Let n G S^+k-i be a paired coloured permutation satisfying 
the hypotheses of the theorem. Starting from n we construct a new sets of 
permutations as follows: 





(b) n F r 



(11) 



S 9 . 



T. 
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1. add the noncoloured element b+ k in any of the b+k possible positions, 
so to obtain b + k new permutations belonging to S^+k] 



2. add the two coloured elements b + k and b + k + 1 in any possible 
positions: this can be done in (b + k)(b + k + 1) different ways, and 
produces permutations belonging to SVffc+l- 

Moreover, the additional hypothesis that the subpermutation of n con- 
stituted by the elements 1, 2, . . . , b — 1 must be isomorphic to r implies that 
r is the minimal permutation of our class, and any other permutation ob- 
tained using the above described construction must avoid the pattern r in 
its elements 1, 2, . . . , b — 1. It is now easy to recognize that, if a permutation 
7r G SVffc-i is given label (b + k), the above construction can be described 
by the mixed succession rule fi given in (I10p . which is enough to conclude. 



Corollary 5.3 7/6=1, then f n = fn = n\F n , and the sequence (f n )neN 
enumerates the class of paired coloured permutations. 

Proof. Set b = 1 in the previous proposition. ■ 

This last sequence also appears in [Slj (it is essentially sequence 
A005442), and can be obtained as the row sums of a particular convolu- 
tion matrix (see [K]). The combinatorial interpretation of (/ n )neN reported 
in [Slj seems to be essentially different from the one given here: it would be 
interesting to have a bijective argument explaining how to relate these two 
interpretations. 

6 Final remarks 

In the present paper we have studied doubled mixed succession rules, 
and, in the commutative case, we have been able to give an expression for 
the sequence associated with one of such rules in terms of the sequences 
associated with the constituent simple succession rules. The next step should 
be to have an analogous result for more general kinds of doubled mixed 
succession rules. For instance, one could consider two succession rules whose 
associated rule operators obey some weaker form of commutativity, such as 
LM = qML, for a given scalar q (or, more generally, LM = f(q)ML, for 
some polynomial /). 

Another presumably fertile line of research concerns exhaustive gener- 
ation. Similarly to what has been done for classical succession rules, one 
can try to develop general exhaustive generation algorithms based on mixed 
succession rules, maybe finding a new way of defining general Gray codes 
depending only on the form of the mixed succession rule under consideration. 
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